Observational Tests of Modified Gravity 
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Modifications of generai reiativity provide an alternative explanation to dark energy for the ob- 
served acceleration of the universe. Modified gravity theories have richer observational consequences 
for large-scale structure than conventional dark energy models, in that different observables are not 
described by a single growth factor even in the linear regime. We examine the relationships be- 
tween perturbations in the metric potentials, density and velocity fields, and discuss strategies for 
measuring them using gravitational lensing, galaxy cluster abundances, galaxy clustering/dynamics 
and the ISW effect. We show how a broad class of gravity theories can be tested by combining 
these probes. A robust way to interpret observations is by constraining two key functions: the ratio 
of the two metric potentials, and the ratio of the Gravitational "constant" in the Poisson equation 
to Newton's constant. We also discuss quasilinear effects that carry signatures of gravity, such as 
through induced three-point correlations. 

Clustering of dark energy can mimic features of modified gravity theories and thus confuse the 
search for distinct signatures of such theories. It can produce pressure perturbations and anisotropic 
stresses, which breaks the equality between the two metric potentials even in general relativity. 
With these two extra degrees of freedom, can a clustered dark energy model mimic modified gravity 
models in all observational tests? We show with specific examples that observational constraints 
on both the metric potentials and density perturbations can in principle distinguish modifications 
of gravity from dark energy models. We compare our result with other recent studies that have 
slightly different assumptions (and apparently contradictory conclusions). 

PACS numbers: 98.65.Dx,95.36.+x,04.50.+h 



I. INTRODUCTION 

The energy contents of the universe pose an interesting 
puzzle, in that general relativity (GR) plus the Standard 
Model of particle physics can only account for about 4% 
of the energy density inferred from observations. By in- 
troducing dark matter and dark energy, which account 
for the remaining 96% of the total energy budget of the 
universe, cosmologists have been able to account for a 
wide range of observations, from the overall expansion of 
the universe to the large scale structure of the early and 
late universe 

The dark matter/dark energy scenario assumes the va- 
lidity of GR at galactic and cosmological scales and in- 
troduces exotic components of matter and energy to ac- 
count for observations. Since GR has not been tested 
independently on these scales, a natural alternative is 
that the failures of GR plus the Standard Model of par- 
ticle physics imply a failure of GR. This possibility, that 
modifications in GR at galactic and cosmological scales 
can replace dark matter and/or dark energy, has become 
an area of active research in recent years. 

Attempts have been made to modify GR at galactic [2| 
or cosmological scales 0, 0, HI- Modified Newtonian 
Dynamics (MOND) and its relativistic version (Tensor- 
Vector- Scalar, TeVeS) are able to replace dark matter 
at galaxy scales to reproduce the galaxy rotation curves, 
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which provided the earliest and most direct evidences for 
the existence of dark matter. The DGP model in 
which gravity lives in a 5D brane world, naturally leads 
to late time acceleration of the universe. Adding a cor- 
rection term f(R) to the Einstein-Hilbert action [|| also 
allows late time acceleration of the universe to be real- 
ized. 

In this paper we will focus on modified gravity (MG) 
theories that are designed as an alternative to dark en- 
ergy to produce the present day acceleration of the uni- 
verse. In these models, such as DGP and f(R) mod- 
els, gravity at late cosmic times and on large-scales de- 
parts from the predictions of GR. We will consider the 
prospects of distinguishing MG models containing dark 
matter but no dark energy from GR models with dark 
matter and dark energy. By design, successful MG 
models will be indistinguishable from viable DE mod- 
els against observations of the expansion history of the 
universe. To break this degeneracy, observations of large- 
scale structure (LSS) must be used to test the growth of 
perturbations. 

LSS in MG theories can be more complicated to pre- 
dict, but is also richer because different observables like 
lensing and galaxy clustering probe independent per- 
turbed variables. This differs from conventional DE sce- 
narios where the linear growth factor of the density field 
fixes all observables on sufficiently large-scales. One of 
the goals of this study is to examine carefully what vari- 
ous LSS observables measure once the assumption of GR 
(with smooth DE) is dropped. 

Structure formation in modified gravity in general dif- 
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fers @, 0, S S M HO, El, El El M M El O from that 
in GR. Theories of LSS in these modified gravity models 
are still in their infancy. However, perturbative calcula- 
tions at large scales have shown that it is promising to 
connect predictions in these theories with observations 
of LSS. Most studies have focused on probes of a single 
growth factor with one or a few observables. In this paper 
we will consider a variety of LSS observables that can be 
measured with high precision with current or planned 
surveys. Our emphasis will be on model-independent 
constraints of MG enabled by combining different ob- 
servables. 

Carrying out robust tests of MG in practice is chal- 
lenging as in the absence of a fundamental theory, the 
modifications to gravity are often parameterized by free 
functions, to be fine tuned and fixed by observations. 
Given the parameter space available to both DE and 
MG theories, it is unclear how the two classes of theories 
can be distinguished. Kunz and Sapone [19( presented 
a rather pessimistic example. They found that one can 
tune a clustered dark energy model to reproduce obser- 
vations of gravitational lensing and matter fluctuations 
in the DGP model. It is not clear if this conclusion ap- 
plies to all modified gravity models and if adding more 
LSS observables helps to break this severe degeneracy. 

In this paper, we first discuss ways of parameterizing 
modified gravity models and dark energy models. §11 
presents the definitions and evolution equations for per- 
turbations in the metric and the energy momentum ten- 
sor. We then classify independent LSS observables based 
on the perturbations that are probed by them. §111 is 
devoted to the use of observational probes of LSS for 
testing MG. We consider the four fundamental pertur- 
bation variables and the observations that can be used 
to probe them. The additional information available in 
the quasilinear regime is discussed in the Appendix. In 
§IV we consider the question of distinguishing MG from 
DE scenarios. The specific question we want to answer 
is: given a set of LSS observations, can a general MG 
model be mimicked by a DE model? If not, what LSS 
observables are required to break the degeneracy? We 
conclude in §V. 



II. PERTURBATION FORMALISM 

By definition, the dark sector (dark matter and dark 
energy) can only be inferred from their gravitational con- 
sequence. In general relativity, gravity is determined 
by the total stress-energy tensor of all matter and en- 
ergy {G^ v = 8wG T^v). Thus we can treat dark mat- 
ter and dark energy as a single entity, without loss of 
physical generality [H [2l], |22 |. This entity has total 
mean matter density pgr and equation of state parame- 
ter w — Pgr/ Pgr- However, when discussing perturba- 
tions in this entity, we may separate it into a matter com- 
ponent (dissipationless particles which can be described 
as a pressure-less fluid free of anisotropic stress) and a 



dark energy component. Throughout this paper, when 
we refer to "smooth" or "clustered" dark energy, we re- 
fer to this dark energy subset of the overall dark sector. 

We may consider the Hubble parameter H(z) to be 
fixed by observations. In a dark energy model, pgr is 
given by the Friedman equation of GR: pgr = 3H 2 /8irG. 
The equation of state parameter is w = — 1 — 2H/3H 2 . 

The corresponding modified gravity model has mat- 
ter density pmg to be determined from its Friedman-like 
equation. We will consider MG models dominated by 
dark matter and baryons at late times and denote fluid 
variables such as the density with subscript mg- 



A. Metric and fluid perturbations 

With the smooth variables fixed, we will consider per- 
turbations as a way of testing the models. In the Newto- 
nian gauge, scalar perturbations to the metric are fully 
specified by two scalar potentials ip and <f>: 



ds 2 = -(1 + 2ip) dt 2 + (1 - 2(j>) a 2 (t) dx 2 



(1) 



where a{i) is the expansion scale factor. This form for 
the perturbed metric is fully general for any metric the- 
ory of gravity, aside from having excluded vector and 
tensor perturbations (see [24[ and references therein for 
justifications). Note that ip corresponds to the Newto- 
nian potential for the acceleration of particles, and that 
in General Relativity <fi — tjj in the absence of anisotropic 
stresses. 

A metric theory of gravity relates the two potentials 
above to the perturbed energy-momentum tensor. We in- 
troduce variables to characterize the density and velocity 
perturbations for a fluid, which we will use to describe 
matter and dark energy (we will also consider pressure 
and anisotropic stress below). The density fluctuation S 
is given by 



6(x,t) 



= P(x,t) - p(t) 



(2) 



where p(x, t) is the density and p(t) is the cosmic mean 
density. The second fluid variable is the divergence of the 
peculiar velocity 



V j T->/(p + p)=V-v, 



(3) 



where v is the (proper) peculiar velocity. Choosing in- 
stead of the vector v implies that we have assumed v to 
be irrotational. This approximation is sufficiently accu- 
rate in the linear regime, even for unconventional dark 
energy models and minimally coupled modified gravity 
models. 

In principle, observations of large-scale structure can 
directly measure the four perturbed variables introduced 
above: the two scalar potentials ip and 0, and the den- 
sity and velocity perturbations specified by S and 8. As 
shown below, these variables are the key to distinguishing 
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modified gravity models from dark energy. Each has a 
scale and redshift dependence, so it is worth noting which 
variables and at what scale and redshift are probed by 
different observations. It is convenient to work with the 
Fourier transforms, such as: 



6{k,t) 



d 3 x 8(x, t) e 



-ik-x 



(4) 



When we refer to length scale A, it corresponds to a 
a statistic such as the power spectrum on wavenumber 
k = 2n/\. We will henceforth work exclusively with the 
Fourier space quantities and drop the " symbol for conve- 
nience. 



B. Evolution and constraint equations 

We consider here the fluid equations for DE and MG 
scenarios. We work in the Newtonian gauge and follow 
the formalism and notation of [201 ] . except that we use 
physical time t instead of conformal time. We are inter- 
ested in the evolution of perturbations after decoupling, 
so we will neglect radiation and neutrinos as sources of 
perturbations. We will make the approximation of non- 
relativistic motions and restrict ourselves to sub-horizon 
length scales. One can also self-consistently neglect time 
derivatives of the metric potentials in comparison to spa- 
tial gradients. These approximations will be referred 
to as the quasi-static, Newtonian regime. We will not 
consider the evolution of perturbations on super-horizon 
length scales; [23j show that differences in their evolution 
may have observable consequences for some MG models 
(discussed further under the CMB below). 



We have allowed for anisotropic stress sources in the en- 
ergy momentum tensor, parameterized by the scalar tr, 
which enters the Euler equation. 

Note that the above equations describe the multi- 
component fluid of baryons, dark matter and dark en- 
ergy; the density and velocity variables for this fluid are 
subscripted gr above (these variables will represent a 
fluid with no dark energy for MG theories below). The 
metric potential variables are <fi and tp in either case. Fur- 
ther, we do not subscript Sp and a as these sources occur 
only in the DE plus GR scenario. 

The linearized constraint equation gives the Poisson 
equation for weak field gravity: 



k 2 <j) = — 47rGa 2 pGR 



<5 G r + 3(1 + w)Ha 



~ -4nGa 2 pGR S G r. 



&GR 
k 2 



(7) 



where in the second line we have dropped the H9GR/k 2 
term as it is negligible for nonrclativistic motions on 
scales well below the horizon. 

Non-zero anisotropic stress a leads to inequality be- 
tween the two potentials: 



ip) = 127rGa 2 (l + w)p a 



(8) 



It is common to take (j> = ip for ordinary matter and 
dark matter; however clustered dark energy can have a 
non-negligible anisotropic stress. 

Eqns. OH] fully describe the evolution of perturbations 
in DE scenarios in the quasi-static, Newtonian regime. 
Next we consider the analogous relations for modified 
gravity scenarios. 



Dark Energy with GR scenario 



2. Modified Gravity scenario 



We first consider the DE scenario, assuming GR. Us- 
ing the perturbed field equations of GR to first order 
gives a set of constraint and evolution equations. The 
evolution of the density and velocity perturbations in- 
cludes gravity and pressure perturbations Sp as sources. 
In the Newtonian limit they give the familiar continuity 
and Euler equations for a perfect fluid. Keeping all first 
order terms, and using the notation S = dS/dt, gives: 



Sgr 



-(l + w)( 
-(l + w)- 



a 
'gr 
a 



30) - 3H — + 3Hw6 GR 



3H 



Sp 



P 

3HwS G r . 



(5) 



In the second line we have dropped the <j> term as it is 
negligible compared to the other terms in the quasi-static 
regime. The Euler equation is given by 



Ogr = -H{1-3w)9gr- 



-0 G r- 



Sp/p 



a 
(6) 



For minimally coupled gravity models with baryons 
and cold dark matter, but without dark energy, we can 
neglect pressure and anisotropic stress terms in the evo- 
lution equations to get the continuity equation: 



'MG 



'MG 



30 



'MG 



(9) 



where the second equality follows from the quasi-static 
approximation as for GR. The Euler equation is: 



Omg 



-HQ 



k 2 yb 



MG 



(10) 



For a generic MG theory, the analog of the constraint 
equations and © can take different forms. We will 
attempt to characterize the general behavior in the weak 
field limit for small perturbations (small S) and non- 
relativistic motions. On sub-horizon scales the field equa- 
tions in MG theories can then be significantly simplified. 
We parameterize modifications in gravity by two func- 
tions G c g(k, t) and j](k, t) to get the analog of the Poisson 
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equation and a second equation connecting <p and ip [251 ] . 
We first write the generalization of the Poisson equation 
in terms of an effective gravitational constant G e s: 

k 2 cj> = ~ATTG c s{k,t)p MG a 2 5MG ■ (H) 

Note that the potential <\> in the Poisson equation comes 
from the spatial part of the metric, whereas it is the 
"Newtonian" potential ip that appears in the Euler equa- 
tion (it is called the Newtonian potential as its gradient 
gives the acceleration of material particles). Thus in MG, 
one cannot directly use the Poisson equation to eliminate 
the potential in the Euler equation. A more useful ver- 
sion of the Poisson equation would relate the sum of the 
potentials, which determine lensing, with the mass den- 
sity. We therefore introduce G e g and write the constraint 
equations for MG as 

k 2 (i, + ^) = -8vrG cff (k, t)p MG a 2 S MG (12) 

<t> = 1>n(k,t) (13) 

where G c g — G e g(l + r)~ 1 )/2. Note that if one starts in 
real space then the corresponding parameters would not 
be Fourier transforms of r\ and G c ff. Thus the Fourier 
transform of the PPN parameter 7 = (f)/ip, the ratio of 
the metric potentials in real space constrained by solar 
system tests, is given by a convolution of rj and ip [8lj |. 
Only if 77 is scale independent would it be the Fourier 
transform of 7. A similar reasoning applies to G c g in 
using the Poisson equation. We prefer to work in Fourier 
space because of the ease of describing perturbations: 
each Fourier mode evolves independently in the large- 
scale, linear regime. Furthermore, the equations describ- 
ing cosmological perturbations in MG theories such as 
f(R) gravity and DGP are generally expressed in Fourier 
space. 

The parameter G c ff characterizes deviations in the 
(ip-\-<p)-5 relation from that in GR. Since the combination 
ip+(p is directly responsible for gravitational lensing, G c g 
has a specific physical meaning: it determines the power 
of matter inhomogeneities to distort light. This is the 
reason we prefer it over working with more direct gener- 
alization of Newton's constant, G e ff. 

The Goff-77 parameterization is equivalent to the Q-t] 
parameterization independently proposed by [18J, where 
Q parameterizes deviations in Poisson equation IJ7J from 
GR. For minimally coupled gravity models, with no dark 
energy fluctuations, it is also equivalent to that proposed 
by [161 ] . And ij is also equivalent to the parameter w pro- 
posed by [13] • DGP and f(R) gravity can be described by 
our parameterization. So as the widely adopted Yukawa 
potential. An exception to our approach is TeVeS as it 
includes scalar and vector fields that are coupled to the 
growth of scalar perturbations. 

For a generic metric theory of MG, one would expect 
that a Poisson-like equation is valid to leading order in 
the potentials and the density perturbation, at least on 



large scales in the linear regime where Fourier modes are 
uncoupled. In this regime, we expect that since the left- 
hand side of the field equations involve curvature, it must 
have second derivatives of the metric perturbations, while 
the right hand side is simply given by the energy mo- 
mentum tensor. On smaller scales, in general a MG the- 
ory may not obey superposition and require higher order 
terms and higher derivatives of the potentials. Similarly 
a generic relation between <p and ip is likely to have a lin- 
earized relation of the form in Eqn. [T21 While it is not 
necessary that the leading term be linear in both the po- 
tentials, observational constraints require that it be very 
close to linear with r\ ~ 1 on small scales where tests of 
gravity exist (see [26[ for a review). 

With the linearized equations above, the evolution of 
either the density or velocity perturbations can be de- 
scribed by a single second order differential equation. In 
the case of MG theories, this equation is simpler as the 
only source is provided by the Newtonian potential ip. 
From Eqns. O and [TU] we get, for the linear solution, 

S(k, t) ~ S 'initial (k)D(k, t), 

$ + 2H5+—t-=0. (14) 
a 

For a given theory, Eqns. [12] and [T3] then allow us to 
substitute for tp in terms of S to determine D{k 1 t), the 
linear growth factor for the density: 

D + 2HD- ^Lp MG a 2 D = 0. (15) 

0- + v) 

We can also use the relations given above to obtain the 
linear growth factors for 9 and the potentials from D. 
Note that in general the growth factors for the potentials 
have a different k dependence than D. In the Appendix 
we give details on the linear and second order solutions 
and summarize quasilinear signatures of MG theories. 

C. Power spectra 

Before we turn to large-scale structure observables, we 
define the power spectra of the perturbed variables. The 
three-dimensional power spectrum of 5(k, z) for instance 
is defined as 

(S(k,z)6(k',z)) = {2vf5v{k + k!)P s {k,z). (16) 

where we have switched the time variable to redshift z. 
The power spectra of perturbations in other quantities 
are defined analogously. We will denote the cross-spectra 
of two different variables with appropriate subscripts, for 
example P$ip denotes the cross-spectrum of the density 5 
and the potential ip. 

We write down next the relation between the power 
spectra of the two potentials and the density in DE and 
MG scenarios. From the Poisson equation ([7]) for GR we 
have 

GR : P (fc, z) = (^G) 2 a 4 p 2 GR Ps > G l{ k ' z) , (17) 
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where P§ is the power spectrum of the potential <p. Using 
the Friedman equation for GR the above equation is often 
written as 



GR:P,(k,z) = iW^WM) 



i 2 k 4 



(18) 



where Hq is the present day value of the Hubble param- 
eter, and Q is the dimensionless density parameter. 

The Poisson equation (fT2)) for MG gives the following 
equations for the power spectra of the metric potentials. 



MG : P^(k,z) = [^G cS {k,z)} 2 a 4 p 2 MG PsMG{k,z)/k A 

(19) 



P m 



[%irG cS {k,z)] 2 4 _ 2 Pg,MG(k,z) 
Y2 a Pmg~ 



[l + rrHKz)}' 



k 4 



where we have used Eqn. rj3]to get the equation for P^. 

For LSS observables, we will need to the power spec- 
tra of (ip + 4>) for lensing, of ip for dynamics, and of S 
for tracers of LSS. We will use Eqns. [T7lfl9l above to 
connect them, along with the relations between the two 
potentials (Eqn. [8] for GR and Eqn. [13] for MG). With 
these relations we can express different observable power 
spectra in terms of a single density power spectrum - for 
MG this will involve the functions G c g(k, z) and rj{k, z). 



III. LARGE-SCALE STRUCTURE 
OBSERVATIONS 

We will assume that the background expansion rate is 
determined by a set of observations: Type la supernovae, 
baryon acoustic oscillation (BAO) and other probes at 
low redshift and the CMB and nucleosynthesis at high 
redshift. These observations measure the luminosity or 
angular diameter distance at a given redshift. The dis- 
tance measures in a spatially flat universe are, within 
factors of 1 + z, simply the comoving coordinate distance: 



X(z) = 



dz' 



(20) 



Furthermore, BAO can directly measure the Hubble con- 
stant at the redshift of galaxies. 

We are interested in the constraints available on per- 
turbed quantities. Hence we will consider observational 
probes of large-scale structure to constrain modified grav- 
ity scenarios. In nearly all cases we will be interested in 
scales in the range 1 — 10 3 Mpc. The MG theories of in- 
terest must modify gravity on horizon scales of order 10 
Mpc; it is an open question how they transition to GR 
on very small scales to satisfy experimental constraints 
from solar system tests. We will assume that the MG 
theories of interest differ from GR over the observation- 
ally accessible scales. 

The most stringent current tests of gravity come from 
laboratory and solar system tests and from binary pulsar 
observations - see [261 ] for a review. Interesting probes 
of gravity on sub-Mpc scales also exist: galaxy rotation 



curves, satellite dynamics, strong lensing observations of 
galaxies and clusters, and X-ray plus lensing observations 
of clusters (e.g. |13|). Modifications in gravity can affect 
the propagation of gravitational wave. Future gravita- 
tional wave experiments such as LISA can detect gravi- 
tational wave from distant supermassive black hole pairs 
in the coalescence phase and thus test this effect [28| . 
We will not consider these tests in this paper. We will 
restrict our attention to large-scale structure on scales 
where theoretical predictions can be made using linear 
or quasilinear perturbation theory. 



Connection of observables to perturbation 
variables 



In principle, observations of large-scale structure can 
directly measure four fundamental variables that describe 
the perturbed metric and (fluid) energy-momentum ten- 
sor: the two scalar potentials i\) and <fi that characterize 
the metric, and the density and velocity perturbations 
specified by 5 and 9. Next we discuss the prospects for 
different probes of these variables. 

Sum of potentials i/j + <f>: Gravitational lensing in 
either the weak or strong lensing regime probes the sum 
of the metric potentials. We will consider the weak lens- 
ing shear (or equivalently the lensing convergence) power 
spectrum as the primary statistical discriminator of MG 
via lensing. 

The spatial components of the geodesic equation for 
a photon trajectory x M (A) (where A parameterizes the 
path) is: 



d\ 2 pa dX dX 



(21) 



For the metric of Eqn. [IJ this gives the following relation 
for the first order perturbation to the photon trajectory 
(generalizing for example from Eqn. 7.72 of [3ll|): 



d 2 x^ 
dX 2 



= -q 2 V^ + (/>) . 



(22) 



where q is the norm of the tangent vector of the unper- 
turbed path. This gives the deflection angle formula 



(23) 



where s = qX is the path length and a.i is the i— th com- 
ponent of the deflection angle (a two-component vector 
on the sky). Since all lensing observables are obtained by 
taking derivatives of the deflection angle, they necessar- 
ily depend only on the combination if> + (j) (to first order 
in the potentials). 

For weak lensing tomography we use the shear power 
spectrum for two sets of source galaxies with redshift 
distributions centered at Zi and Zj. Following standard 
treatments of weak lensing, this may be derived from the 
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deflection angle formula to get the shear power spectrum 
on angular wavenumber I (|32[): 

<W0 = fdxW i (x)W j (x)k- 4 P i ,+Jk=-,x) , (24) 



where the weight function Wi is simply 



Xi - X 
Xi 



(25) 



for source galaxies at a single comoving distance Xi = 
x( z i) (it can be easily generalized for sources specified by 
a redshift distribution). We have assumed a flat back- 
ground geometry for simplicity; our results throughout 
this paper can be generalized to a curved spatial geome- 
try by replacing x m the argument of W by the angular 
diameter distance. 

Note that in the literature the lensing power spectra 
for GR are expressed in terms of the density power spec- 
trum Ps(k) assuming the standard Poisson and Friedman 
equations. Usually anisotropic stress is neglected so that 
one can substitute into the above equation the relation 
between the power spectra: = 9 k~ 4 HQ Q 2 P$/a 2 

from Eqn. [TSJ For MG, this substitution breaks down 
due to the modifications of the Poisson equation and the 
Friedman equation. However the correct substitution can 
be made in terms of G e s(k, z) using equation [T9l and the 
modified Friedman equation (which depends on the spe- 
cific theory). 

Since lensing probes the sum of the metric poten- 
tials, with the deflection angle formula following from 
the geodesic equation (which simply describes how cur- 
vature affects trajectories), it may not by itself test the 
field equations of the gravity theory. However lensing 
measurements at multiple source redshifts are sensitive 
to the growth of the lensing potential, which does of- 
fer a test of the MG theory. And by combining lensing 
with other observables, the relation of P^,+<j, to P$ can be 
tested. Recent studies that have examined constraints on 
MG theories with weak lensing include [IMl2lll3^. l34l . l35l ]. 

Another important observable in lensing is galaxy- 
galaxy lensing, the mean tangential shear around fore- 
ground (lens) galaxies. Its Fourier transform, the galaxy- 
lensing cross-spectrum, depends on tp + <j) and on the 
galaxy number density. It is given by an equation similar 
to Eqn. 1241 with the power spectrum of the lensing po- 
tential in the integrand replaced by the three-dimensional 
cross-power spectrum, and with one of the weight func- 
tions replaced by one representing the foreground galaxy 
distribution: 



dx W ^^ P g( , +d k=i, X ), (26) 

"-A X 



where W g i is the normalized (foreground) galaxy redshift 
distribution (e.g. [36j]). Galaxy-galaxy lensing has been 
well measured from the SDSS survey. It is a very useful 
check on galaxy bias, hence it aids the interpretation of 
galaxy clustering measurements ((37[) as well. 



Assumptions: In using weak lensing observations with 
the above formalism, one assumes that intrinsic corre- 
lations are negligible or removable (in general these can 
differ for different gravity theories), that the weak lensing 
approximation is valid, and that galaxy properties that 
affect photometric redshift determination are not affected 
by the gravity theory. 

Newtonian Potential ip: This can be measured by 
dynamical probes, typically involving galaxy or cluster 
velocity measurements. If gravity is the only force deter- 
mining galaxy accelerations at large scales (as expected) , 
we have from Eqn. 1101 



k 2 tp = 



d(a9 g ) 
dt 



(27) 



where 6 g = V-v 9 . On sub-Mpc scales this relation can be 
used to constrain ip using galaxy satellite dynamics and 
rotation curves (e.g. [SH). Redshift distortion effects in 
the galaxy power spectrum probe larger scales, which we 
address in more detail here. 

The redshift space power spectrum of galaxies is a well 
measured quantity. It can be expressed in the large-scale, 
small angle limit as (e.g. [3!|) : 



W) = 



2u 2 u* 
P a {k) + ^ r P a e g {k) + —Pe g {k) 



F 



H 2 {z) 

where u = kn/k is the cosine of the angle of the k vector 
with respect to radial direction; P g , P g g , Pe g are the real 
space galaxy power spectra of galaxies, galaxy-6* ff and 9 g , 
respectively; a v is the ID velocity dispersion; and F(x) 
is a smoothing function, normalized to unity at x = 0, 
determined by the velocity probability distribution. The 
dependence on u enables separate measurements of all 
three power spectra, though Pg is the hardest to measure 
with high precision [40l |41| . Furthermore, measurements 
of P g at smaller scales provide information on pairwise 
velocity dispersion a v [43| . 

In the linear regime, we can rewrite Eqn. 1271 as 



fc 2 V> 



d{aD e )/dt^ 
Dg ^ 



(29) 



Here Dg is the growth factor of 9 g . For MG models, 
Dg has a simple relation to D, the linear density growth 
factor: Dg oc aD = a(3HD 1 where (3 = din D/dlna. 
In the linear regime we have 9 g (k, t) = 9 g (k.,ti)Dg(k.,t). 
Note that the above equation does not require Dg to be 
scale independent, so it is applicable to modified gravity 
models and clustered dark energy models. Note also that 
we do not distinguish the growth factor of 9 g from that of 
9 because we only use its time (redshift) derivative, which 
is expected to be very similar. Velocity measurements at 
multiple redshifts are required to measure ip from the 
above equation, as described in 42 1. 

For clustered DE models, the galaxy v g is not neces- 
sarily equal to v of the total fluid. From the Euler equa- 
tion ([(j]) applied separately to different components of the 
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fluid, we can see that the DM and DE velocities evolve 
differently since only the latter is affected by pressure 
perturbations in the DE. As a first order approximation, 
galaxies and baryonic gas velocities trace that of the DM. 
So what one actually measures is 9 g ~ 9dm ^ $de ^ 0. 
This distinction can be relevant for DE models with large 
perturbations on sub-horizon scales if these are not corre- 
lated with the matter fluctuations (i.e. if the DE power 
spectrum has a different shape from the matter power 
spectrum) . 

Assumptions/Caveats: The galaxy peculiar velocity 
only probes ip where there are galaxies. So potentially 
there is a bias related to the environment of galaxies. 
However, since gravity is a long range force, the poten- 
tial where galaxies reside is determined by matter over 
a much larger region and thus should be unbiased with 
respect to the overall ip. Galaxies themselves are not 
sufficiently massive to contribute to this long range po- 
tential. However, to obtain v g from limited redshift bins, 
one does need to parameterize the redshift dependence 

of Vg. 

The accuracy of the velocity information inferred from 
the redshift space galaxy power spectrum relies on the 
modeling of the redshift distortion. The derivation of 
Eqn. [28] is quite general - it can be applied to general DE 
or MG models. However, Eqn. [55] does not describe red- 
shift distortions to percent level accuracy 39]. Nonethe- 
less, with improved modeling of the correlation function 
in redshift space [H| the associated systematic errors in 
velocity (and tp) measurements can be reduced. 

Density contrast 5: The clustering of galaxies is one 
of the earliest measures of large-scale structure, and its 
measurements have advanced over the last three decades. 
The galaxy power spectrum P g is the simplest statistical 
measure of correlations in the galaxy number density. 
Several other probes of large-scale structure also probe 
the density field: clustering of the Lyman-alpha forest, 
clustering of quasars and galaxy clusters, the abundance 
of galaxy clusters, and (in the future) 21-cm emission 
measurements of the high-redshift universe. 

However, given a measured galaxy power spectrum, 
the power spectrum P$ of the underlying mass density 5 
may differ due to galaxy bias. Further the galaxy-density 
relation may be non-local and vary slightly in different 
gravity theories due to differences in the tidal field that 
influence collapsed objects such as galaxy halos. We will 
restrict ourselves to large scales (fc k a \, the nonlinear 
wavenumber) where bias is scale independent in simple 
models of galaxy formation. This allows us to infer the 
mass power spectrum from the galaxy power spectrum 
without detailed modeling of their relation, because it 
is possible to fit for the bias directly from the data. We 
discuss below the caveats to this assumption for clustered 
dark energy. 

The galaxy density in three-dimensional space may be 
expressed in terms of the density and bias parameters 61 



and 62 as 

5 g = *± = bl 8 + %P. (30) 

Tig Z 

This expansion is useful for small values of 5; it can be 
used in a perturbative expansion to explore what mea- 
surements are sufficient to measure the bias parameters 
61 , bi as well as 5 (see [45| for details on the bias formal- 
ism) . Eqn. [2FJ above shows how the three-dimensional 
galaxy power spectrum P g can be obtained from redshift 
space measurements. A second way of measuring P g is 
from imaging data with photometric redshifts. This pro- 
vides measurements of the angular power spectrum of 
galaxies, which is a projection of the three-dimensional 
galaxy power spectrum 

r W 2 (y) ( I \ 

Cg(l)= Jd X ^P g (k=-,x), (31) 

where W g is the normalized redshift distribution of galax- 
ies included in the sample. With good photo-z's it is a 
narrow range with width of order 0.1 in redshift, so that 
many such angular spectra can be measured at different 
mean redshifts from a survey (e.g. [Ifjj]). 

1. Galaxy bias with clustered dark energy 

In clustered dark energy models it is not a priori clear 
whether the galaxy overdensity is related to the matter 
overdensity 8 m or to the total fluid overdensity <5gr- We 
argue below that at least for some galaxy populations, 
S g is directly related to 5 mi even though the evolution 
of the matter density responds to the full gravitational 
potential (which receives contributions from dark energy 
clustering as well). 

One way to see this is to consider the centers of mass 
of galaxy halos at sufficiently high redshift Zi that the 
dark energy density is negligible. The clustering of these 
halo centers is then simply a biased version of the mass 
distribution. Hence at z\ one can write S g (zi) = b(zi)5M, 
with b(zi) independent of scale for large enough scales. 
As they evolve to redshifts below unity, their motions 
are given by the potential if), just as for the matter 
field. Hence their evolution obeys the continuity and Eu- 
ler equations: S g ~ ~9 g /a and 9 g ~ —H9 g + k 2 ip/a. 
The matter density obeys the same equations with 8m 
and 9m as the density and velocity perturbations. This 
means that the bias factor preserves its scale indepen- 
dence: at low redshift, it relates the galaxy power spec- 
trum to the matter power spectrum and is not directly 
sensitive to the clustering of dark energy. For example 
the halo model expression [4?], [H[ for the bias evolu- 
tion is: b(z) ~ 1 + (y — 1)/S sc (z) where S sc (z) cx D(z) 
is the density required for spherical collapse at 2, and 
v = S sc (z)/a with a the smoothed rms mass fluctuation. 
The expression for ellipsoidal collapse has two additional 
parameters but still has no scale dependence. 
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Clustered dark energy follow Eq. [6] with w ^ and 
a 7^ 0, so it has a different time and spatial dependence 
from 8 m . If the dark energy clusters significantly, it is 
therefore possible that galaxies have a scale dependent 
bias relative to it and therefore to the total density field. 

The above argument is very general but relies on some 
approximations. These are well justified for massive ha- 
los, for which the evolution at low redshift is very sim- 
ple: consider galaxy halos of mass M 3> M», where M* 
is the standard halo model nonlinear mass. The centers 
of mass of these halos can be mapped to high-cr peaks in 
the nearly Gaussian mass distribution at high redshift. 
Moreover, they do not move significantly, so it is evident 
that their power spectrum at large scales evolves simply 
by the growth of its amplitude. Such massive halos corre- 
spond to galaxy clusters and LRG's at moderate to high 
redshift. For galaxies in lower mass halos, halo motions 
and mergers change their clustering at low redshift, so 
one has to be careful in modeling their bias factors. 

Another route to 6 in any GR scenarios is through the 
metric potentials. Given lensing measurements of tp + <f> 
and dynamical measurements of ip, one can obtain the 
potential <j>. Using this, the Poisson equation then 
gives 5, since the Gravitational constant is known in GR. 
Thus S is not independent of the metric potentials even 
for clustered DE models. 



2. Empirical determination of bias 

To leading order then, knowledge of b\ allows us to 
relate P g to P$. Barring extreme scenarios of clustered 
dark energy, we take 5 to be the full density field. 

Provided a halo-model description applies reasonably 
well to our universe, bias can be determined by combin- 
ing observations and using two and three-point statistics. 
For concreteness we consider the bias parameters b\ , 62 
that can be determined from the data using the power 
spectrum and bispectrum (denoted B) measurements. In 
a deterministic bias model, one can then get the density 
power spectrum. With P g = b\P$ and the reduced three- 
point parameter Q ~ B/P 2 (see the Appendix and [45[ 
for full expressions), one has a relationship between the 
Q parameter of galaxies and mass 



By using P g and measurements of Q g for different tri- 
angles, both bias parameters and Ps can be determined. 
(A similar analysis can be done in real space, e.g. using 
counts in cells. The skewness S3 is given by the shape 
of the power spectrum and bias parameters.) While this 
is a simplified model, it helps us address what changes 
for MG: the predictions for Pg and Qs both change, with 
the former given by the new linear growth factor on large 
scales and the latter by next order terms in perturbation 
theory (see the Appendix for more details) . For well spec- 



ified gravity scenarios, these calculations can be done and 
thus the bias factors determined from measurements. 

A second approach to measuring b\ is to use the galaxy- 
mass cross-correlation measured by galaxy-galaxy lens- 
ing in combination with the galaxy power spectrum (e.g. 
[37I]). This has the advantage that one uses only two- 
point statistics that can be measured with high accuracy. 
However, as discussed below and by [25[, for MG theo- 
ries there is a complication because the Poisson equation 
is needed as well since lensing measures the potentials 
rather than 5. So for MG theories, the extraction of the 
bias parameter in this approach is more complicated - 
but nevertheless feasible by jointly fitting for bias and 

G e ff- 

3. Galaxy cluster mass function 

A different probe of 5 is provided by the mass func- 
tion of galaxy clusters. Given Gaussian initial conditions 
and a spherical/ellipsoidal collapse model, the number 
density of galaxy clusters can be related to the linear 
density contrast. In the spherical collapse scenario, a re- 
gion containing mass M will collapse if the overall density 
fluctuation exceeds a threshold S c . The number of such 
regions can be predicted from the Gaussian statistics and 
this fixes the halo mass function dn/dM, the number of 
halos with mass M. 

In the standard ACDM cosmology, gravitational dy- 
namics is determined by GR. The mass function of galaxy 
clusters is sensitive to the smoothed mass density vari- 
ance a\ on scale R, which is dependent on the cluster 
mass and is typically of order 10 Mpc (e.g. [HI]). This is 
related to the density power spectrum as: 

/d 3 k 
^3 P s (k)W t 2 op _ hat (kR), (33) 

where Wt op -hat is the window function for averaging with 
a spherical top-hat. 

For clustered DE models, the cluster formation picture 
becomes complicated. The presence of the anisotropic 
stress invalidates the spherical collapse model and more 
complicated models such as ellipsoidal collapse with tidal 
fields need to be used. Furthermore, the fate of an over- 
dense region is no longer determined by the matter fluc- 
tuation 5 m alone. DE fluctuations <5de and a affect %j) 
through equations [7] and [3 And Sp affects the evolution 
of Sue through Eqns. [S]and[ni Thus a combination of S m , 
<^DE: Sp and a act in determining the evolution of a given 
region of matter - the resulting collapse condition has 
yet to be worked out. Since many galaxy clusters form 
recently at z < 1, where DE is non-negligible, DE fluctu- 
ations could leave some detectable signatures in cluster 
abundance. 

For probing the dark universe, this is a valuable fea- 
ture. It implies that galaxy cluster abundances contain 
information on not only fluctuations in matter but also 
fluctuations in dark energy, and thus is a promising probe 
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of the total 5. To extract such information requires fur- 
ther modeling of the DE model, but a simplified model 
can be obtained as follows. The collapse condition based 
on energy conservation should be linear in the DM and 
DE perturbation variables, since they are all first order 
variables of the energy-momentum tensor. At high red- 
shift, the dark energy contribution should vanish (assum- 
ing pde pm)- Thus we may assume that matter fluc- 
tuations are the only source of growth for the late time 
S m as well as <5de, dp and o~ responsible for the LSS. In 
this picture, all perturbation variables are correlated and 
have deterministic relations [52]. The collapse condition 
can be simplified into a modified condition on S m alone. 
An effective 5° can be defined for specific DE models, 
such that when a region reaches S m > , it will collapse. 

The usual collapse model deals with isolated objects 
and thus BirkhofFs theorem is implicitly required. Mod- 
ifications in GR result in a generic breakdown of the 
Birkhoff's theorem. This significantly complicates the 
modeling of cluster abundance in MG models, since the 
fate of a given region is determined not only by matter 
and energy inside this region, but also matter and energy 
outside. However, given a MG model, one can still pre- 
dict the probability for a given region with overdensity 
<5 TO to collapse and thus predict cluster abundances. 

Assumption/Caveats: Unlike the use of gravitational 
lcnsing to probe ij) + it is model dependent to probe 
5 from cluster abundance. (1) The cluster abundance 
requires careful modeling, even in the simplest case of 
smooth DE models. For example, the tidal field makes 
the spherical collapse model only a rough approximation. 
(2) The observable-mass relation is needed to connect 
observable (e.g. X-ray flux, SZ flux or cluster richness) 
to the mass of clusters. These cluster properties often 
involve complicated gastrophysical processes and can not 
be predicted with sufficiently high precision from first 
principles. As a consequence, using cluster abundance to 
probe S often require model-dependent calibrations. 

In spite of these caveats, it may be hoped that the well 
posed problem of the evolution of a region in an initially 
Gaussian random field will be calculable, and related to 
the linear density field in generic MG or DE models. 

Velocity divergence 9: Many existing velocity mea- 
surement are based on distance indicators: the difference 
of the true distance from what is inferred from the re- 
cession velocity gives an estimate of the peculiar velocity 
of a sample of galaxies or clusters [53]. The pairwise 
velocity at small separation can be measured through 
anisotropic galaxy clustering in redshift space at cosmo- 
logical distances [431] . While challenging, there are ongo- 
ing attempts to improve measurements of bulk flow mea- 
surements, based on SNe la [54|. An independent method 
is the kinetic Sunyaev Zel'dovich (SZ) effect [55[ of clus- 
ters which is directly proportional to the cluster peculiar 
velocity and enables a rather model independent mea- 
surement method [561 ] . These measurement are likely to 
have lower signal-to-noise than the redshift space distor- 



tions discussed above. Further it is unclear whether they 
estimate 6 of the total fluid in a clustered DE scenario, 
for the reason discussed above. 

CMB: The CMB power spectrum is given by: 

Ctt(1) = JdkJ d x ' F CMB (k, I, x ') ji[k X (z')] (34) 

where the spherical Bessel function ji is the geometric 
term through which the CMB power spectrum depends 
on the distance to the last scattering surface. The func- 
tion Fcmb combines several terms describing the pri- 
mordial power spectrum and the growth of the potential. 
We will regard Fcmb as identical to the GR prediction 
since we do not invoke MG in the early universe (up to 
the redshift at last scattering). 

The CMB anisotropy does receive contributions at red- 
shifts below last scattering, in particular due to the in- 
tegrated Sachs Wolfe (ISW) effect [13] • In the presence 
of dark energy or due to modifications in gravity, gravi- 
tational potentials are in general time varying and thus 
produce a net change in the energy of CMB photons: 

AT f djijj + <j>) 

— Iisw = - / J t a Kt)dx • (35) 

The ISW effect, like gravitational lensing, depends on 
and probes the combination if; + <fi. The ISW signal is 
overwhelmed by the primary CMB at all scales (although 
it does produce a bump at the largest scales in the CMB 
power spectrum). For this reason, it has to be measured 
indirectly, through cross-correlation with other tracers 
of large scale structure. The resulting cross-correlation 
signal is then 

<WQ = /^ + 0)( fc = ^)« 2 ^- 06) 

Here, P g fj )+ ^\ (k,\) is the cross-power spectrum of (ip + 

<f>) and galaxies or other tracers of the LSS such as 
quasars or clusters. By cross-correlating the CMB tem- 
perature with galaxy over-density S g , the ISW effect has 
been detected at < 5a confidence level [58| and pro- 
vides independent evidence for dark energy, given the 
prior of a spatially flat universe and GR. This cross cor- 
relation signal depends on galaxy bias, which has to be 
marginalized to infer cosmology. With the aid of gravita- 
tional lensing, uncertainties of galaxy bias can be avoided 
[HIH^I- Furthermore, since the ISW amplitude peaks on 
the largest scales, it also has a strong correlation with 
large scale bulk flows and produces a cross correlation 
signal with potentially better signal-to-noise than that 
of the density-ISW cross correlation [63] . 

The primary CMB is Gaussian and statistically 
isotropic. However, gravitational lensing distorts the 
CMB sky and induces anisotropy and Fourier mode- 
coupling in the CMB, which should not exist otherwise. 
This feature should allows reconstruction of the lensing 
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potential from future high resolution CMB maps [64j . 
The CMB sky is the furthest lensing source and thus can 
probe tfj + (f> at redshifts well above unity. This will be 
useful to constrain those MG and DE models in which 
deviations from ACDM persist at these redshifts. 

ISW measurements and future measurements of lens- 
ing and galaxy clustering can probe scales approaching 
the horizon scale. This provides an additional test of 
MG models in which the growth of perturbations is al- 
tered at relatively high redshift on super-horizon scales. 
[2! ] showed that growth on super-horizon scales is con- 
strained to be universal for MG models with ip = <p. [23| 
show how it differs for f(R) models which do not obey 
this constraint, and describe the transition from super- 
horizon to sub-horizon scales. If measurements achieve 
high accuracy on these large scales, they can be combined 
with information on sub-horizon scales to provide addi- 
tional constraints on the ratio of potentials 77 for such 
MG models. 

Summary: The quantity that can be measured most 
robustly is the sum of potentials ip + <f>, through gravita- 
tional lensing and the ISW effect. With a bit more mod- 
eling, the Newtonian potential tp can be inferred from 
galaxy velocity measurements (i.e. redshift space distor- 
tions). To obtain model independent constraints on the 
total density perturbation S is challenging if one allows 
for dark energy clustering in the GR scenario. Galaxy 
clustering is likely to be an effective measure of the mat- 
ter fluctuation S m , while cluster abundance is a promis- 
ing probe of 6 as it is sensitive to DE fluctuations as well. 
Although the galaxy peculiar velocity is likely to be well 
measured in the future, the DE peculiar velocity (and 
therefore the overall v and 6) is likely the most difficult 
to measure. Cross-correlations of large-scale structure 
tracers with the lensing potential or the Newtonian po- 
tential are probably the most promising tests of MG in 
the near future, as we discuss next. 



B. Joint constraints from multiple observations 

If multiple observables are to be combined, model in- 
dependent information can only be inferred if they probe 
the same range of redshift and length scale. The distance- 
redshift relation will be measured to ~ 1% accuracy by 
the next generation SNIa and BAO surveys at low-z and 
by the CMB at high-z. The next generation BAO sur- 
veys can further measure H(z) at low redshift. With 
the expansion rate of DE and MG models tightly con- 
strained, measurements of perturbed variables become 
powerful discriminators. 

The distance-redshift relation at redshift z is given by 
an integral over the expansion rate, and therefore the en- 
ergy densities, from redshift to z. This measurement 
at z £ 1 has provided evidence of acceleration, consistent 
with ACDM. On the other hand, CMB measurements at 
high-z for both distances and perturbations are consis- 
tent with a universe governed by GR, with its energy 
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FIG. 1: Upper panel: We plot the normalized distance d(z) 
(solid curves, almost coincident) and the linear growth rate 
D(z)/a(z) (dotted curves) for < 2 < 3 for two dark energy 
models (in black and red) and a DGP model (in blue). The 
distances and growth rates are normalized to give 1 at high 
redshift (z=1100). Lower panel: The fractional deviations 
in distance (solid curves) and linear growth (dotted curves) 
from the fiducial ACDM model are shown for a dark energy 
model (black) and DGP (blue) (see text for details). Note 
that the DGP growth curve is sensitive to the parameters 
chosen to fit the distance redshift relation: for a distance 
curve that matches ACDM better, the growth curve would 
also have smaller deviation. 

density dominated by matter and radiation (36|. Thus 
either dark energy or modification of gravity must pro- 
duce effects that are significant at z 1 and negligible at 
z ~ 1000. In Fig. [I] we show as examples the deviation 
(from ACDM with tt de = 0.7) of a model with tt de = 0.75 
and a flat DGP model with il m = 0.3. It is clear that for 
both distances and perturbations, significant deviations 
occur at low-z in such models [82| . 

The most promising scale/redshift range in the near 
future is ~ 10— 100s Mpc at redshifts ~ 0.3 — 1. Imaging 
and spectroscopic observations are likely to be made on 
these scales and will be robust to many sources of error 
and dependence on specific models. We list below sev- 
eral categories of surveys that will test MG and DE mod- 
els. Two sets of surveys are indicated: surveys planned 
for the near future (significant data within 5 years), and 
surveys planned to start in about a decade. (The list is 
not complete as several projects have been formulated or 
modified recently.) 

• Multicolor imaging survey: With photometric red- 
shifts for millions of galaxies, these surveys pro- 



11 



vide measurements of weak lensing, galaxy cluster 
abundances, and the angular clustering of galaxies, 
clusters and quasars. These measurements probe 
ip + 4>, S g and their cross-correlation. Upcoming 
surveys include: 

- DES, KIDS, PS1, HSC (2008-). < z < 1. 

- LSST, SNAP, DUNE (2014-). < z < 3. 

• Spectroscopic surveys: While primarily designed to 
measure the distance-redshift relation and H(z) us- 
ing the baryon acoustic oscillations in galaxy power 
spectra, they will provide improved measurements 
of P g , P gv , P vv on large scales. Some surveys will 
target z ,$ 1 galaxies and others will select galaxies 
at higher rcdshift, 2 ^ z ,$ 3. 

- LAMOST, WiggleZ, HETDEX, WFMOS, 
BOSS (2008-). < z < 3. 

- ADEPT (2014-). 1 < z < 2. 

• 21 cm surveys: SKA [H[(2015-). The square kilo- 
meter array (SKA) has the potential to detect ~ 1 
billion galaxies over < z £ 1.5, with a deeper sur- 
vey extending to z ~ 5, through 21cm line emis- 
sion of neutral hydrogen in galaxies. If success- 
ful, it will provide high precision measurements of 
the distance-redshift relation through BAO's 65], 
and tests of MG through: (a) weak lensing maps 
with accuracy comparable to that of large optical 
surveys (66[, (b) velocity measurements through 
redshift distortions of galaxy clustering [25], and 
(c) ISW measurements through CMB-galaxy cross- 
correlations. 

• SZ and X-Ray Cluster surveys: These will measure 
the abundances of galaxy clusters out to z ~ 1 
and beyond. Cosmological applications will depend 
on supplementary optical data to get photometric 
redshifts of the detected clusters. 

- SZA, SPT, ACT, APEX, eROSITA (2008-). 
< z < 1. 

• CMB: temperature and polarization maps provide 
high-z constraints and also measurements of the 
ISW effect and CMB lensing, which are probes of 
if) + (f> at lower redshift. 

- PLANCK and ground based missions (2008-) 

We have indicated the approximate redshift range over 
which these surveys will provide accurate measurements. 
It would be most useful to have different observables over- 
lap in redshift and length scale in the range z ~ 0.3 — 1 
and at scale A ~ 10 to several 100 Mpc. This range of 
scales covers the linear and quasilinear regimes of struc- 
ture formation (we are assuming that MG effects are 
present on these scales). We consider next two promis- 
ing combinations of observables that on the 5-7 year 



timescale will enable measurements of the MG functions 
G e ff and r\ on these scales. 

Lensing and galaxy power spectra: Planned next- 
generation imaging surveys (see above) will have area 
coverage in excess of 1000 sq. degrees, enabling few per- 
cent level measurements of lensing power spectra. The 
same imaging surveys will also measure the angular clus- 
tering of the galaxies C g (at z ~ 0.3 — 0.6, the redshift 
of the lensing mass) to percent level accuracy; cluster 
abundances will also be measured through optical and 
SZ surveys: both measurements probe the matter den- 
sity S. Alternatively, spectroscopic surveys like BOSS 
will measure P g , the three-dimensional power spectrum, 
to percent level accuracy. The shear power spectra can 
be combined with the density power spectra measured at 
z ~ 0.3 — 0.6 and scales of 10-100 Mpc. Using the Poisson 
equation, G e g will then be tightly constrained, assuming 
statistical errors dominate the error budget. The main 
galaxy sample and LRG sample of the SDSS has already 
been used in constraining MG models through galaxy 
clustering alone (e.g. (6?| ), though in a model-dependent 
way. 

Cross-correlations of galaxies with shear and 
velocity: Galaxy-galaxy lensing measurements made 
from imaging surveys probe the lensing potential-galaxy 
correlation. This measurement has been made to high 
accuracy from the SDSS [H, [6i|. In the near future one 
can expect measurements of P gv to a few percent (e.g. 
from the BOSS survey) at z ~ 0.3 — 0.6. In combination 
with percent-level galaxy-galaxy lensing measurements 
over the same range of redshift, the ratio of potentials r\ 
will be precisely constrained [421 ] . 

The measurements described above would be major 
advances in constraining MG theories, as the current con- 
straints on ~ 10 — 100 Mpc scales are weak and insuffi- 
cient to test MG theories with any robustness. For par- 
ticular models the scale and redshift evolution of a single 
statistic, such as the lensing power spectrum, can be pow- 
erful as well. We leave for future work a detailed study 
of how well these measurements will test MG theories. 
In considering an observable suitable for distinguishing 
models of gravity, one must address the familiar prob- 
lems in extracting cosmological information due to sta- 
tistical and systematic errors, i.e. the expected precision 
on the measurement, the physical assumptions necessary 
to connect observable to the four variables of interest, 
and the degeneracy with other cosmological parameters. 



IV. MODIFIED GRAVITY VS. DARK ENERGY 

Specific models of MG and DE can be tested by com- 
bining observations of the expansion rate and large scale 
structure. For example, in the ACDM model (and well 
defined scalar field models), the growth of the large scale 
structure is completely determined by the expansion his- 
tory: there exists a fixed relation between the expansion 
rate and the growth of LSS. This consistency check has 
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been carried out in the literature and is indeed able to dis- 
tinguish specific models investigated Furthermore, 
this consistency check can be performed in a model inde- 
pendent way to search for signatures of violation of GR, 
with the prior of smooth DE [70] . Current data pass this 
consistency check, although violations of the consistency 
relation can not be ruled out [7l[ . 

DE models that depart from scalar field models can be 
much more complicated, with a break down of the cor- 
respondence between the expansion rate and large scale 
structure. The expansion rate is determined by pgr and 
w. However, two extra DE properties, the anisotropic 
stress a and the response of pressure to perturbations 
(Sp), can affect the growth of the LSS. These two proper- 
ties are determined by the microphysics of the DE model 
and are independent of pgr and w. As a consequence, the 
growth of LSS is no longer fixed by the expansion rate 
and the above consistency check can not be applied to 
search for signatures of the violation of GR. Current ob- 



servational constraints on the sound speed 
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[72| and the anisotropic stress [73j are weak. Further- 
more, these studies use a particular form of a and 6p 
and assume that one can be switched off when studying 
the other. Thus a potentially wide range of DE models 
with non-negligible anisotropic stress and pressure fluctu- 
ations are still viable against observations. To investigate 
the feasibility of distinguishing between DE and MG, we 
will allow for arbitrary anisotropic stress and pressure 
perturbations. 

Modifications of gravity (at least the class of theories 
we have considered) involve two extra quantities which 
govern LSS, namely, modification of Newton's constant, 
G g, and the ratio of potentials 77 = <j>/ip. Although 
these quantities determine the gravitational interaction 
of perturbations, they do in general affect the expansion 
rate H(z) - unlike for GR and its Newtonian limit. This 
is in part due to the fact that for MG Birkhoff 's theorem 
no longer holds and thus the usual exercise of calculating 
H(z) from the Newtonian dynamics of a spherical matter 
distribution no longer applies. Thus G c ff and 77 represent 
real extra degrees of freedom in MG theories. 

The extra degrees of freedom in MG and clustered DE 
models can produce similar observational consequences. 
For example, the anisotropic stress breaks the equality 
between <p and ip, mimicking the role of 77 in MG mod- 
els. Thus one might expect that by tuning the two extra 
degrees of freedom in DE models, one can mimic a given 
MG model to fit observations. Indeed, Kunz & Sapone 
[l9l ] explicitly construct a DE model which reproduces de- 
generate^, ip and S m with the flat DGP modified gravity 
model. [84J 

This degeneracy certainly deserves further investiga- 
tions. In this section, we consider in more details the 
question: can one always succeed in tuning DE models to 
produce observational consequences identical to a given 
MG model? If the answer is yes, then one can never 
unambiguously test for deviations from GR. 

The answer to the above question is incomplete in 



fully describing the dark degeneracy. The complemen- 
tary question, which needs to be answered is: can one 
always tune MG models to produce observational con- 
sequence identical to a given DE model? If the answer 
is yes, then one can never unambiguously justify the ex- 
istence of DE. However, this question is more difficult 
since it requires a general parameterization of the rela- 
tion between the expansion rate H(z) and the nature of 
a general gravity theory - such a parameterization is not 
yet available. This limit in theoretical understanding of 
MG forces us to investigate only the first question, since 
we know the most general way of parameterizing the in- 
fluence of DE on the expansion history of the universe. 
Furthermore, we have constrained our study to a special 
class of MG models, in which gravity is minimally cou- 
pled to matter. The study of both questions for the most 
general MG models is beyond the scope of this paper. 

The relationship between the four perturbation vari- 
ables (j>, if), S and 9 is fixed for a complete DE or MG 
theory. These consistency relations are the key to prob- 
ing the nature of DE and MG. With just two variables 
being observable, one can only test against one consis- 
tency relation and, as we see below, by tuning the two 
extra degrees of freedom in clustered DE models, any 
MG model can be mimicked. However, with more ob- 
served variables, one can test other consistency relations 
and hope to break the degeneracy between DE and MG 
models. In this section, we explore the feasibility of dis- 
tinguishing DE and MG models. In this section we con- 
sider only the question of distinguishability in principle, 
without regard to the accuracy of observations in the 
foreseeable future. 



A. Two perturbation observables 

First we assume that both potentials are observables, 
i.e. we require <f> and ip to be identical in the two mod- 
els. From the discussions in previous sections, these two 
quantities are the most likely to be measured to high pre- 
cision. So we set them identical in the constraint equa- 
tions for GR and MG to get relations between the re- 
maining variables. Comparing Eqn. [8] for a in GR with 
the constraint Eqns. [T2l and IT51 for MG gives 



2rj- 
a = — - 



1 G c ff Pmg 



3 r\ 



1 G pgr 



<5mg 



(37) 



In addition by combining the Poisson equation |(7J) for 
GR with Eqns. [12] and [13] for MG, we obtain a second 
constraint 



S GR + 3(1 + w)hJ^ = 



if 



2 G c ff pmg 
1 + 1 G "pgr"' 



JMG 



(38) 



The question then is whether Eqns. [5] [6] and [38] have 
solutions for 5gr, $GR and Sp, in terms of MG variables 
(recall that a is now fixed by Eqn. |3"T|) . Without a 
fundamental theory, Sp can take any form (2l| : hence 
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there is always a form of Sp satisfying all three equa- 
tions. Namely, there is always a DE model which can 
mimic the given MG model to produce identical <j> and 

v>. 

The degeneracy persists for other combinations of two 
perturbation variables. We have discussed above that in 
a clustered DE model, it is difficult to establish whether 
galaxies, other tracers or cluster abundances probe <5 or 
5 rn (or neither!). If we assume that a sub-set of LSS 
observations will provide measurements of 5, then com- 
bining that with measurements oitp + <p from lensing, we 
have 



GW pmg 
3(1 + to)" I G pgr 



-S 



(39) 



Thus if a is free, it can be chosen to match the above 
equation for any set of theories. 

Extra information can break this degeneracy. The 
response of pressure to density perturbations and the 
anisotropic stress are determined by the microphysics 
of the DE model. It requires a theory to provide such 
closure relations (see [23| for more detailed discussions). 
For example, the quintessence model predicts vanishing 
a and negligible pressure perturbation on sub-horizon 
scale. Even if advances in the understanding of general 
DE theory do not provide such specific information, some 
general constraints can still break the degeneracy. For ex- 



ample, if Sp takes the form Sp = c 2 Sp and c 2 



;{t), as 



is true for the adiabatic case, solutions do not exist in 
general for equations [5] [H and [3S] In this case, one can 
not find a DE model to mimic the given MG model. 

Another physically well motived example is for the 
anisotropic stress a. A natural source of a is the velocity 
perturbations in the fluid. By the requirement on gauge 
invariance, the evolution in a may be parameterized in 
the following form in the Newtonian gauge [lU, [71] , 



a + 3H& = 



3 1 



(40) 



where c V j S is the viscous parameter. This equation in 
general contradicts equation [37J and [35] above and thus 
no DE model that satisfied Eq. [40] can mimic the given 
MG model. 

Extra information can also come from additional ob- 
servables. The equations above show that if we have just 
one additional observable, such as S or 0, there will in 
general be no solution for the remaining two variables 
that satisfies three equations (e.g. [5] [6] and [38]). We 
consider this next. 



B. Three or more observables 

If both potentials and S are observable then the theory 
is constrained much more tightly, especially if they are 
measured multiple redshifts. 
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FIG. 2: First consistency condition for at least one DE model 
to mimic <j>, ip and S in a flat DGP model. The dashed line 
given by Eqn. [41] represents the required condition, while the 
solid curve is the actual relation in flat DGP. When a — > 0, 
77 — > 1 and for a — > 00, 77 — > 1/2. The points on the curve 
with a — 0.5 and a = 1 are indicated. (For flat DGP lines 
with different £l m lie on top of each other.) The disagreement 
between the two curves shows that DGP is a modified gravity 
model that can not be mimicked by any dark energy model. 



For a DE model to mimic the given MG model, S, </> 
and ip must satisfy the three equations [7J [T^] and Q2] 
This imposes the the following consistency relation for 
G e ff and 77, 



V 



, G c ff pmg 
■~G 



1 



Pgr 



(41) 



So if the given MG model does not obey the above rela- 
tion, no DE model can produce S, <j) and ip identical to 
the given MG model, no matter how the DE properties 
are fine-tuned. 

Eq. 0J] represents a strong constraint on MG models 
as it shrinks the 2-parameter rj-G c s space in MG models 
into a straight line. We show as an example that the 
DGP model does not satisfy this condition. 

In a flat DGP model, Goff = G and 77" 1 = (1 + 
1/3#dgp)/(1-1/3/3dgp) EH- Here/W = l-2r c JT(l+ 
H/3H 2 ) = 1 - 2r c H/{\ - 2r c H) < 0, with H 2 = 
H/r c + fl m a~ 3 and r c = 1/(1 — f2 m ). We have normal- 
ized H(z = 0) = 1. For the DGP model, p M G = Q m a~ 3 
(up to a normalization, which is irrelevant for this dis- 
cussion). By requiring the DE model to reproduce the 
expansion history of the given MG model, we have Pgr = 
H 2 . Fig. [2] shows that the consistency condition is sig- 
nificantly violated by all flat DGP models. This means 
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Combining both constraints on <r s , we obtain 

tT 1 ~w/3 + 1 + 3w . (46) 

w is fixed by the condition pgr = H 2 . (3 is calculated 
from the given MG theory. So the above equation can 
be checked from the viewpoint of MG models unambigu- 
ously 

Again, Fig. [3] shows that this condition is severely vio- 
lated for the DGP model: thus no DE model can mimic 
a flat DGP model to reproduce identical <fi, ip and 5. We 
have verified that this is true for f(R) models as well. 

These relations present general constraints, without re- 
sort to real observation data. Observations show that at 
the present epoch, w < —1/3 since the universe is ac- 
celerating, while (3 > since structure is growing. From 
Eq. 251 we have c 2 < —1/3, if the related DE model can 
reproduce cf>, ip and 8. Furthermore, Eq. 25] tells that 
1] < today. 

V. DISCUSSION 



FIG. 3: The second necessary condition for at least one DE 
model to mimic (f>, tp and 5 in the given DGP model. The 
condition r/" 1 ~ w(3 + 1 + 3w, given by Eqn. 1461 implies the 
variable on the y-axis should be zero for all a. For flat DGP 
with Q. m = 0.2, this condition is also severely violated for 
a > 0. 



that no DE model can produce 0, ip and 5 identical to 
a flat DGP model that satisfies observational constraints 
on the expansion history. 

This conclusion seems to contradict [l{|. However, [l9| 
require the dark matter density fluctuation S m to be iden- 
tical in the GR and MG scenarios. We require the total 
6 to be identical instead. In GR, the Poisson equation 
specifies the cj)-S relation, not the (j)-S m relation, so with 
S as an observable, one directly construct consistency re- 
lations and thus distinguish between DE and MG. 

There is also a constraint on the anisotropic stress, 
from Eq. 1 [Ml HI and [H 



V 



- 1 



8 . 



3(1 + w) 

Comparing Eq. [5] with Eq. I10[ we obtain 



9 ZHw 



w 



1 



1 



) - 

10 la 



Since HO ~ l3aH 2 8 < k 2 S/a, we have 

2 (1 + w)<7 r\~ x — 1 
° s ~ 5 = 3 
Comparing Eq. Owith Eq. [51 we obtain 



2 w(6/a-Z$) (3 
°^ -3H5 3+1) 



(42) 



(43) 



(44) 



(45) 



We have described the role of perturbations in test- 
ing theories of modified gravity (MG) against large-scale 
structure observations. We have chosen the class of MG 
theories that are described by scalar perturbations, so 
that two metric potentials suffice to describe the per- 
turbed space-time. We then consider the quasi-static, 
Newtonian limit of the perturbed field equations and 
compare dark energy (DE) and MG theories. 

Our main focus is on the relationship of different 
observables - lensing, large-scale dynamics of galaxies, 
galaxy clustering, cluster abundances and various cross- 
correlations - to the four perturbation variables of MG 
theories: the two metric potentials, the density field, and 
the divergence of the peculiar velocity. In mil] we give 
the relationship of measured power spectra in real space, 
redshift space and on the sky with theoretical predictions 
for these perturbation variables. We highlight the use of 
two effective functions to test MG theories: the ratio of 
the metric potentials and the effective Gravitational 
constant G c g. We also consider in the Appendix quasi- 
linear signatures of MG and show how the MG functions 
affect second order corrections to the power spectrum 
and the bispectrum. 

We discuss in detail what is actually measured by var- 
ious large-scale structure observations once the assump- 
tions of smooth dark energy and GR are dropped. While 
lensing and dynamical probes have a direct connection to 
different potential variables, tracers of the density field 
and cluster abundances must be treated carefully in ex- 
tracting information about the density field from them. 
The most robust tests of MG effects can be made by 
combining different observables from planned multi-color 
imaging surveys and redshift surveys: see ElIIIBI for two 
examples that will be feasible in the near future. 

Observables that may not be useful in constraining 
smooth DE can be crucial for testing MG because there 
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are more variables to be measured and different observ- 
ables are sensitive to them. For instance, the redshift 
space power spectrum P g g is not considered a valuable 
probe of structure formation in DE studies as other meth- 
ods produce lower statistical errors on dark energy pa- 
rameters. But in a MG scenario, P g g is useful because 
it probes the Newtonian potential tp that other probes 
are not sensitive to. It can be combined with galaxy- 
galaxy lensing, which probes P g t^, +( j,\, to constrain the 
ratio of potentials (p/ip f §III Bj) . This is a case where 
observables from multi-color imaging surveys and spec- 
troscopic surveys must be combined to test MG theories. 
More generally, once one allows for MG scenarios, mul- 
tiple observables are needed to test theories. Thus the 
diversity of LSS observations, which has lost some of its 
appeal in the recent trend of going after a single dark en- 
ergy figure of merit, becomes vital (see [59| for a broader 
criticism of dark energy driven research, and 60] for a 
rebuttal). 

Finally we consider a question posed recently in the 
literature: can a DE model be constructed to mimic any 
MG theory? We show that with observations of multi- 
ple perturbed variables (three are sufficient in general), 
unique signatures of MG theories can be established. 
We show with the example of the flat DGP model how, 
given sufficiently accurate measurements of the lensing 
and Newtonian potentials and the density, no DE model 
can mimic the DGP model. 

Our results may be compared with other recent stud- 
ies in the literature, some of which appeared while this 
work was in progress [H, HE [Z^l ■ These studies tackle 
the question of how dark energy and modified gravity can 
be distinguished. We have clarified the appar ent conflicts 
between this paper and [H in ffVBl 0,1751, [ll present 
a formal argument that any deviation from GR can be 
absorbed into the dark sector as an effective dark com- 
ponent. The effective stress-energy tensor of this compo- 
nent, T°„, is defined as the deviation of the given MG 
theory from Einstein's field equations. Tj^ is conserved, 
as expected for the usual DE model. From this argument, 
one might conclude that MG can not be distinguished 
from DE gravitationally. However, [23| also pointed out 
that the effective dark component has a generic, though 
implicit, coupling to matter, despite the conservation of 
T^. This coupling hides in the closure relations for this 
dark component, which depend on external matter and 
the metric, instead of just its internal microphysics. A 
theory with such a dark component is fundamentally dif- 
ferent from conventional DE models of the kind we have 
considered, even ones with strong clustering of the DE. 
(Indeed, it should not be surprising that allowing for a 
dark component with an arbitrary stress-energy tensor 
and couplings to matter can mimic any modification of 
gravity.) Our result, that DE models - with no coupling 
to matter - can be distinguished from MG, appears to 
be consistent with the analysis of [23]. In i]IVI we showed 
how consistency relations, obtained from the evolution 
and constraint equations obeyed by perturbation vari- 



ables, help to distinguish between DE and MG. The vio- 
lation of the consistency conditions that can occur with 
sufficient observables (see ^IV B[) would thus imply either 
the modification of gravity or a coupling of the "dark en- 
ergy" with matter in a GR scenario. 

There are several assumptions and caveats in this 
study. We study MG theories with scalar perturba- 
tions to the metric; our formalism does not apply to 
theories with additional vector or tensor degrees of free- 
dom. While describing the quasilinear regime of large- 
scale structure (where planned observations have the best 
signal-to-noise), one needs to be aware that it is not clear 
how to obtain nonlinear predictions of some MG theo- 
ries. We have chosen to work in Fourier space, where 
the description of clustering is simpler, but this means 
that there is not always a direct relationship of our ef- 
fective functions for MG theories with the real space de- 
scription of these theories (e.g. G e s is related to its real 
space counterpart in the Poisson equation by a convolu- 
tion with the density field). And finally, we have left for 
future work a detailed study of the accuracy with which 
MG tests can be performed by the next generation of 
surveys. 



VI. APPENDIX: PERTURBATION THEORY IN 
MODIFIED GRAVITY 

The fluid equations in the Newtonian regime are given 
by the continuity, Euler and Poisson equations. Keep- 
ing the nonlinear terms that have been discarded in the 
study of linear perturbations in the rest of the paper, the 
continuity equations gives: 



8 + 9 = - 



(47) 



where the term on the right shows the nonlinear coupling 
of modes. Note that the time derivatives are with respect 
to conformal time in this Appendix. The Euler equation 
is 



d 3 fci fe 2 fci-(fc-fci; 

(2^) 3 2fc2|fc 1 _ jfc2|2 



6 + H6 - k 2 ij) = I _ * a „Zr 0(fci)fl(fc ~ fci) 

(48) 

We neglect pressure and anisotropic stress as the energy 
density is taken to be dominated by non-relativistic mat- 
ter [73]. The Poisson equation is given by Eqn. [T2l and 
supplemented by the relation between ip and given by 
Eqn. 1131 Using these equations we can substitute for %j) 
in the Euler equation to get 



6 + HO 



Pmgo 2 <5 



(l + »7) 
d 3 h fc 2 fci ■ (fc - fci) 
(2tt) 3 2k\\k x - A?2|2 



9(k 1 )9(k-k 1 ) (49) 



Eqns. [17] and 05] are two equations for the two variables 
5 and 9. They constitute a fully nonlinear description 
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of MG theories and can be solved once rj and G e ff are 
specified. An important caveat is that they may never- 
theless be invalid for particular theories, for example if 
the superposition principle is violated. 

Next we consider perturbative expansions for the den- 
sity field and the resulting behavior of the power spec- 
trum and bispectrum. Let S = 81 + 82 + ... where 
62 ~ 0{8\). Higher order effects due to gravitational 
dynamics become detectable on 10s of Mpc at low red- 
shift. While this is strictly true only for general relativity, 
any MG theory that is close enough to GR to fit obser- 
vations can also be expected to have this feature. In the 
quasilinear regime, i.e. on length scales between ^10-100 
Mpc, mode coupling effects can be calculated using per- 
turbation theory. For MG, let us simplify the notation 
by introducing the function: 



equations [47] and [49] as follows 



Cwg(M) 



87rGcff 
(1 + V) ' 



(50) 



which is simply AttG in GR but can vary with time and 
scale in MG theories. The evolution of the linear growth 
factor is given by substituting for ip in Eqn. [51] to get 



Si + HSi - C,MGPMGa 2 Si = 0. 



(51) 



In GR, the relation of ip to S is given by the Poisson equa- 
tion with constant G. In MG, this relation involves both 
•q and G G ff. If either of these functions have a depen- 
dence on k or z, then the solution for the growth factor 
changes. The linear solutions for ip and ip + <f> are then 
simply obtained using Eqns. [T2"l and [TBI 

We show below that in addition the second order solu- 
tion has a functional dependence on G Q g and rj that can 
differ. Thus potentially distinct signatures of the scale 
and time dependence of G e s(k,z) can be inferred from 
higher order terms. These rely either on features in k and 
t in measurements of P$+$ and P$, or on the three-point 
functions, which even at a single redshift can have dis- 
tinct signatures of MG [78[ . Quasilinear signatures due 
to T](k, z) can also be detected via second order terms in 
the redshift distortion relations for the power spectrum 
and bispectrum. Our discussion generalizes that of [fjj 
who examined a Yukawa-like modification of the Newto- 
nian potential. 



A. Second order solution 

From a perturbative treatment of Eqns. H71 and |4"51 the 
second order term for the growth of the density field is 
given by 



S 2 + H5 2 - PMG a ^ Cm g S2 = 
HI 1 [6 1 ,6 1 ] + I 2 [S 1 ,S 1 ]+ ii[5\,Si], 



(52) 



where I\ and I2 denote convolution like integrals of the 
two arguments shown, given by the right-hand side of 



h[Si,Si](k) = 



d 3 ki k ■ ki 



k\ 



6i(ki)Si(k - ki) (53) 



and 



h[SiJi}(k) = 



d 3 k! k 2 h • (jfc-fci) 
(2tt) 3 2kf\ki 



Si(ki)Si(k-h). 

(54) 



Finally, the last term in Eqn. [52] is simply Ii[Si,Si] — 
Ii [Si, Si]+Ii[Si, Si]. Note that by continuing the iteration 
higher order solutions can be obtained. 

From the above equations it follows that if Cmg = 
Cmg(^) then 62 may be specified by /c-integrals over 
Si, so that one may express the functional relationship 
S2 = S2[Si] (where it is understood that 5% at a given 
wavenumber k depends on Si at all other wavenumbers). 
But for the general case of a MG theory with scale de- 
pendent QmGi the second order solution has additional 
scale and time dependence behavior that is not deter- 
mined by the linear solution (owing to the third term on 
the left hand side in Eqn. 152"]) . So the functional relation- 
ship must be modified to: S 2 = 8 2 [5i, Cmg]- This means 
that quasilinear evolution provides an additional signa- 
ture of MG. That is, even if the initial power spectrum 
is not fully specified (e.g. if the running of the spectral 
index is not well constrained), the comparison of linear 
and quasilinear growth rates can reveal the signature of 
MG. In practice, whether the quasilinear signature is sig- 
nificant must be determined by computations for specific 
models (see [7!| for a specific model for which it is not). 
Note also that the second order correction to the density 
power spectrum also involves the third order density field 
as it is given by P2 ~ {S 2 } + (S183). The qualitative fea- 
tures we highlight for 82 will also be found in 83, which 
is also given by iterations of the nonlinear Eqns. [?7] and 
EU 

We summarize the comparison of linear and second 
order solutions for the density for GR versus MG. We 
have identified the function (,MG(k,t) as containing all 
the information about MG that affects density and ve- 
locity fields. For the density field the first and second 
order solutions can be compared to GR as: 

• Linear growth in GR: In smooth dark energy GR 
models, 8i(k,t) is a separable function of scale and 
time. 

• Linear growth in MG: In MG theories, 8i{k,t) is a 
separable function of k and t if and only if the MG 
function Cmg is independent of scale. 

• Second order solution in GR: In smooth dark en- 
ergy GR models, the second order solution S 2 (k,t) 
is not separable. It is however determined by inte- 
grals over 81. 

• Second order solution in MG: In MG models with 
CMc(k,t), 82 is no longer determined solely by Si 
and contains additional signatures of MG. 
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Note that for weak lensing measurements, quasilinear 
corrections are given by the density times G e ff (by substi- 
tuting higher order terms into Eqn[19|). So the resulting 
signatures can be straightforwardly computed using the 
higher order solutions for the density field. 



B. Three-point correlations 

Distinct quasilinear effects are found in three-point cor- 
relations (we will use the Fourier space bispectrum), as it 
is the lowest order probe of gravitationally induced non- 
Gaussianity. The bispectrum for the density field Bg is 
defined by 

(5(ki)6(k2)5(k 3 )) = (27r) 3 fo(fci + + k 3 )B s {h,k 2 , k 3 ). 

(55) 

Since Bg ~ {5 3 ) ~ (<5 2 <5 2 ) (using (5f) = for an ini- 
tially Gaussian density field), the second order solution 
enters at leading order in the bispectrum. Note also that 
the wavevector arguments of the bispectrum form a tri- 
angle due to the Dirac delta function on the right-hand 
side above. In practice, a very useful measure of non- 
Gaussianity is the reduced bispectrum function Q, which 
for the density field 5 is given by 



then has 



Qs = 



B s (ki,k 2 ,k 3 ) 



P6(ki)P S (k 2 ) + Ps(k 2 )P s (k 3 ) + P s (h)P s (k 3 ) 

(56) 

To leading order Q is independent of the amplitude of the 
linear power spectrum (both numerator and denominator 
are 0(<5f), see [78]) and is nearly constant with triangle 
size in GR. It is however sensitive to the shape of the 
triangle. The dependence on size and shape changes for 
MG theories and is in principle a probe of Cmg- It is 
beyond the scope of this paper to elaborate on the mea- 
surement of the bispectrum from galaxy surveys; we will 
instead focus on the prospects for lensing measurements. 
[67l | have tested Yukawa like modifications of gravity us- 
ing Q for the galaxy density measured in real and redshift 
space. 

The lensing bispectrum contains perhaps the clearest 
signature of MG. It is a projection of the three dimen- 
sional bispectrum 

k 6 B^ +4> ~ (8irG cfi a 2 p MG ) 3 {5 3 ) =< (87rG cS a 2 p MG ) 3 (Sfa) 

Since both Si and 6 2 are function of CmGj measurements 
of B^, +l p are sensitive to G e g and Cmg separately. 

The reduced lensing bispectrum in a MG theory can 
be expressed in terms of the density power spectrum and 
bispectrum as: 

Q G c ff(fci)G c ff(fc 2 )G c ff(fc 3 ).B,5(fci,fc 2 , k 3 )/p M G 

^ +<P * fc3 2 G c 2 ff (fc 1 )P 5 (fc 1 )G 2 ff (fc 2 )P 5 (fc 2 )/fc2fc2 +sym ... 

(58) 

For equilateral triangles, Q in MG theories is simpler 
since the G c ft factors in all the terms are the same. One 



„ k Qs(MG) 
Vj/>+</>(MG) -~ Z 

Goft Pmg 



(59) 



The ratio of Q for MG versus GR for equilateral triangles 
is given by 



Qip+4>(MG) Qs(MG) G PGR 
— OC 



b+<j>(GR) 



(GR) G G ff PMG 



(60) 



Note that Qg itself depends on Cmg- Bernardeau [78| 
shows that with n = 1 but a scale dependent G e ff, Qs for 
given initial power spectrum is relatively insensitive to 
Cmg- If that holds for generic intial power spectra and 
gravity models, it would imply that Q^+mmg) probes 
G e ff for models with rj = 1. 

In general a measurement of the lensing power spec- 
trum and reduced bispectrum (roughly speaking, oiP^ + ^ 
and Q^ +( p) is sufficient to measure departures from GR. 
There are three underlying functions (P$, Cmg and G g) 
to be determined. For given k and source redshift, we 
have measurements of P and of Q as a function of tri- 
angle shape. Thus while the equilateral triangles may 
be regarded as sensitive primarily to G e ff, elongated tri- 
angles will be sensitive to Cmg, an d therefore to n. In 
practice one must take account of the fact that the bispec- 
trum has lower signal-to-noise than the power spectrum 
on quasilinear scales (8(j, so one must fit for the desired 
information from all triangle configurations and sizes to 
constrain the MG functions. 

To summarize this section, quasilinear effects thus offer 
two signatures of MG. 

• A scale and time dependent feature on quasilinear 
scales in the power spectrum that depends on -q 
and G e ff. This enters through the second order 
contribution to the power spectrum. 

• Signatures in the bispectrum: additional signatures 
of modified gravity are present in three-point cor- 
relations of the density and potential fields. Inde- 
pendent of the shape and amplitude of the power 
spectrum, the dependence of the reduced bispec- 
trum Q on triangle size and shape is a useful test 
of MG. The reduced lensing bispectrum for exam- 
ple has a strong dependence on G c g. 

We have not considered here whether a clustered DE 
model can mimic both these signatures. It would be of 
interest to carry out the second order calculations for a 
set of MG models and compare predicted deviations with 
observational error bars. 
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